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ABSTRACT 

The origin of ultracompact dwarf (UCD) galaxies, compact extragalactic stellar sys- 
tems, is still a puzzle for present galaxy formation models. We present the compre- 
hensive analysis of high resolution multi-object spectroscopic data for a sample of 24 
Fornax cluster UC Ds obtained with VL T FLAMES. It comprises previously published 
data for 19 objects (jMieske et a l. 2008) which we re-analysed, including 13 with avail- 
able HST photometric data. Using Virtual Observatory technologies we found archival 
HST images for two more UCDs and then determined their structural properties. For 
all objects we derived internal velocity dispersions, stellar population parameters, 
and stellar mass-to-light ratios (M/i)* by fitting individual simple stellar population 
(SSP) synthetic spectra convolved with a Gaussian against the observed spectra using 
the NBuRSTS full spectral fitting technique. For 14 objects we estimated dynamical 
masses suggesting no dark matter (DM) in 12 of them and no more than 40 per cent 
DM mass fraction in the remaining two, in contrast to findings for several UCDs in 
the Virgo cluster. Some Fornax UCDs even have too high values of (M/L)* estimated 
using the Kroupa stellar initial mass function (IMF) resulting in negative formally 
computed DM mass fractions. The objects with too high (M/L)* ratios compared 
to the dynamical ones have relatively short dynamical relaxation timescales, close to 
the Hubble time or below. We therefore suggest that their lower dynamical ratios 
(M/L)^y n are caused by low-mass star depletion due to dynamical evolution. Over- 
all, the observed UCD characteristics suggest at least two formation channels: tidal 
threshing of nucleated dwarf galaxies for massive UCDs (~ 10 s M©), and a classical 
scenario of red globular cluster formation for lower-mass UCDs (<; 10 7 M©). 

Key words: galaxies: dwarf - galaxies: elliptical and lenticular, cD - galaxies: evo- 
lution - galaxies: stellar content - galaxies: kinematics and dynamics 



1 INTRODUCTION 



Ultra-compact d warf galaxies dHilker et al.l 1 19991 : 
iDrinkwater et all l200d : iPhillipps et all 120011 )" initially 
discovered as extragalactic sources unresolved from the 
ground-based observations, represent a new class of com- 
pact stellar systems (CSS) observed in the nearby Universe 



* Based on the archival ESO VLT data (programme 078. B- 
0496) available through http://archivc.eso.org/ and Hubble Space 
Telescope archival data (programmes 8090 and 10129) available 
through http://hla.stsci.edu/ 

f E-mail: igor.chilingarianSSastro.unistra.fr, chil@sai.msu.ru 



JMieske et all 12004 iHasegan et ail 120051; I Jones et all 120061 : 
IMieske et aTp20071 . 12003 : iMisgeld et all I2008T ). At least 
an order of magnitude smaller (10 < R e < 100 pc) than 
M32-like compact elliptical (cE) g alaxies, but still signifi- 
cantly larger than globular clusters l|Drinkwater et a"i1l2003l ; 
I Jordan et ai1l2005l ). UCDs are best studied in the two nearby 



Jordan et al.ll^UUa ). utus are best stu died m tne t wo nearby 
clusters of galaxies : Fornax and Virgo dHasegan et al"] 20051: 



Jones et all 120061: lEvstigneeva et all 120071; iHilker et all 



20071 : IMieske et al.|[200l ). Typically, UCDs have luminosi- 
ties between —13.5 < My < —10.5 mag and masses about 
2 • 10 6 < M < 10 s M (|Mieske et al.ll2008l ). 

UCDs, initially defined on a morphological basis, 
may indeed represent a heterogeneous class of objects 
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(|Mieske et al.l l2006h of different origins. The concepts of 
UCD formation include: (1) very massive globular clus- 
ters h aving the same origin as "normal" ones (jMieske et al.l 
I2002T I, however brightest UCDs (My < -12 mag) are 
too bright and frequent to be statistically accounted for 
by the Gaussian representation of the globular cluster lu- 
minosity function; (2) st ellar superclusters formed ingas- 
rich mergers of galaxies (jFellhauer fc Kroupal I2OO2I , l2005h : 
(3) end-products of small-scale primordial density fi uctu- 
ations in dense environments (|Phillipps et al.l 200 



:t al.1 l200ll ; 



tidally stripped nu cleated dEs (dEN's, iBekki et al.l _ 
iGoerdt et ail |200ej) or simply dE N's with very low surfac e 
brightness outer components (cf. iDrinkwater et al1l200St ). 
None of the concepts is completely ruled out or confirmed 
yet. 

The dynamical M/L ratios of UCDs are on average 
about twice as large as those o f Galactic globular clusters 
of comparable metallicity (e.g. iDabringhausen et al. Il2008l : 
IMieske et aT]|200St iForbes et al.l |200S| ; iTavlor et aljboiol ). 
although environmental differences exist (see below). It has 
been shown that about 1/4 of the M/L offset with respect to 
globular cl usters is probably due to dynamical evolution of 
the latter (jKruiissen fc Mieskej[2009h . The remaining offset 
in M/L indicates that UCDs may mark the on-set of dark 
matter domination in sm all stellar systems (|Gilmore et al.l 
l2007l ; lGoerdt et alj|2008l), o r indeed probe a variation of t he 
IM F (babringhausen etai1l200Sl : IMieske fc Kroupall200sl ) . 

IMieske et al. ( 20081 ) state a need in observations aimed 
at the studies of UCDs stellar populations and age de- 
termination in particular, required to explain the low dy- 
namical M I L ratios in Fornax UCDs compared to those 
in the Virgo cluster. Here we re-analyse the same obser- 
vational data applying a powerful full-spectral fitting tech- 
nique allowing us to obtain simultaneously internal kinemat- 
ics and stellar population properties in the largest UCD sam- 
ple available until now. We also improve some data reduc- 
tion steps, necessary because of the more stringent require- 
ments for stellar population analysis regarding sky subtrac- 
tion. Another motivation for re-analysing these data was to 
double-chec k the velocity disper sion measurements for cer- 
tain UCDs. IMieske et al] (|2008l ) used late-type giant stars 
from uj Cen having metallicities of around —1.0 dex as stellar 
templates. While this is close to the average UCD metallic- 
ity, UCDs span a wider range of metallicities and for metal- 
rich objects the template mismatch due to the metallicity 
difference m ay affect velocity d ispersion measurements as 
suggested in IChilingarianl (|2006h . 

The paper is organised as follows: in Section 2 we de- 
scribe spectroscopic data reduction and analysis, as well as 
complementary archival HST data used to extend the sam- 
ple of UCDs with known structural parameters; the results of 
the full spectral fitting are presented in Section 3; compari- 
son of kinematical and stellar population properties of UCDs 
with early-type galaxies and a general discussion about the 
origin of UCDs is given in Section 4. 



2 DATA: SOURCES, REDUCTION, ANALYSIS 
2.1 Spectral data 

We used the data obtained in the course of the study of 
CSSs in the Fornax cluster (program 078.B-0496, P.I.: L. In- 



fante). The datasets are publicly available through the ESO 
Data Archiv^E The observations were collected at the ESO 
Very Large Telescope with the FLAMES/Giraffe spectro- 
graph (|Pasauini et alj2002h in the multi-object "MEDUSA" 
mode, using the HR9 setup giving a resolving power R m 
17000 in the wavelength range 5120-5450 A in the ser- 
vice mode in 14 observing blocks (OB) between June 
and December 2007 with the total integration time of 
54 000 s ec. All parti c ular details about observations arc 
given m IMieske et al] (|2008f) . We use the distance modulus 
of th e Fornax cluster m — M = 31.39 mag (|Freedman et al] 
l200lf ) corresponding to the distance 19 Mpc and the spatial 
scale 92 pc arcsec - . 

In order to perform the stellar population analysis we 
had to improve the data reduction and sky subtraction with 
respect to the requirements for the kinematical analysis pre- 
sented in IMieske et al] (|2008l ). therefore we had to intro- 
duce several additional steps before and after standard ESO 
FLAMES pipeline data reduction. 

The first additional pre-processing step is to apply 
the L aplacian cosmic ray cleaning algorithm (|van Dokkuml 
120011 ) to the original science frames and mask all the regions 
affected by cosmic ray hits in every individual frame. 

The second step, the diffuse light subtraction, is of a 
great importance for the success of our study. The simulta- 
neous calibration (SimCal) Th-Ar lamp inside the spectro- 
graph was switched on during the observations. The average 
flux level in the SimCal fibres exceeded those of scientific 
targets by at least two orders of magnitude. Therefore, the 
scattered light from bright arc lines at a level of 1-2 per cent 
severely contaminated the neighbouring fibres. There is also 
a smooth diffuse light component in the spectrograph, hav- 
ing however, only little effects on the flat-fielding. 

To account for this contamination, we created a dedi- 
cated scattered light modelling algorithm. Since FLAMES 
fibres are sparsely placed on a CCD plane: the median inter- 
fibre distance is 14.3 pixels while the fibre FWHM is around 
4 pixels, therefore it is possible to use inter-fibre gaps to es- 
timate the scattered light contribution. We use fibre traces 
derived from a flat field to identify inter-fibre gaps which 
normally should not contain any signal. Then, we assume 
that all signal detected in the gaps represent a superposition 
of scattered light and global diffuse light. In the beginning 
we smooth this contribution along dispersion with a boxcar 
of 4 pix, then we fit it across dispersion at every position 
along fibre traces between each pair of SimCal fibres using 
the 3-rd order smoothing splines. We use the spline interpo- 
lation instead of smoothing splines in a small region of the 
frame containing a parasite signal generated by the CCD 
read-out amplifier. 

The two-dimensional scattered light model constructed 
in this fashion is created and subtracted from every spectral 
image. In Fig [1] we present a fragment of a spectral image 
before and after scattered light subtraction. 

Then, the corrected files are used to feed the ESO 
FLAMES data reduction pipeline ESOREX to proceed 
through all steps up-to the linearisation of extracted spectra. 

Then, we apply the third, post-processing step includ- 
ing the adjustment of the wavelength scale to account for 



http: / /archive. eso.org/ 



Masses of Fornax U CDs 3 



J 



50 




200 



Figure 1. A fragment of the science frame before (top) and after 
(bottom) diffuse light subtraction. The dispersion is in the vertical 
direction. The bright trace in the middle is the SimCal fibre, the 
one to the right is one of the brightest scientific targets (F-19). 



heliocentric radial velocity corrections, sky subtraction, and 
combination of spectra for multiple exposures. All further 
data analysis is performed on individual one-dimensional ex- 
tracted spectra of every object. 

We fitted the high-resolution pegase.hr 

|Le Borgne et all l2004h simple stellar population (SSP) 
models against the observational data using the modi- 
fication of the NBursts full spectral fitting technique 
(jChilinearian et all l2007bl lah. The original version of 
NBursts cannot be used directly wit h our data as we 
did fo r FLAMES-LR04 observations (jChilinearian et al.1 
l2008d lbh , because the spectral resolution of FLAMES in the 
HR9 setup exceeds that of the pegase.hr stellar population 
models (R ~ 10000). Therefore, we degraded the spectral 
resolution of the data in order to match that of pegase.hr 
and be able to provide unbiased estimates of stellar velocity 
dispersions. 

We used the high signal-to-noise twilight spectra ob- 
tained in the HR9 setup of FLAMES available for down- 
load from the instrument web-pages. Then, we fitted the 
spectra in 5 segments along the wavelength range in ev- 
ery fibr e against the Solar s pectrum from the elodie.3.1 
library { Prugnicl et al. 2007) using the ppxf method by 
ICappellari fc Emselleml (j2004T ). The procedure was similar 
to the determinati on of the line-spread-functi on (LSF) vari- 
ations described in lChilingarian et al.l (|2008ch used for low- 
resolution FLAMES spectra with the principal difference 
that now FLAMES-HR9 twilight spectra served as "tem- 
plates" and elodie.3.1 spectrum as "data". The derived 
differential spectral line spread is almost constant along the 
wavelength direction and does not exhibit any significant 
fibre-to-fibre variations. Its shape can be described as a 
Gauss-Hermite function with do = km s~\ a — 12 km s _1 , 
/i3 = 0, hi = —0.05. Then, the data convolved with this 
parametrization matched the resolution of pegase.hr mod- 
els. This procedure is very similar to the degradation of the 
pegase.hr resolution we performed in all previous studies 



using the NBursts technique to take into account spectro- 
graph's LSF. 

The contamination of UCD spectra by the light of 
NGC 1399 is an important source of biases for the stel- 
lar population parameter estimates. As it creates an ad- 
ditive background with very smooth spectral features due 
to the high intrinsic velocity dispersion of the galaxy (a > 
200 km s _1 ), the metallicity measurements can become un- 
derestimated while velocity dispersion measurements be- 
come o verestimated (see Appendix A3 in lChilingarian et al.l 
l2007bl . for details). These effects should become very impor- 
tant in the inner region of the cluster. The sky fibres were 
placed far away from the galaxy centre and were hardly a 
subject to any contamination by the NGC 1399 halo, there- 
fore we do not expect any over-subtraction for UCDs located 
at large projected distances from NGC 1399. 

To account for NGC 1399's halo contamination, 
we used its spectrum at the position (a,<5)./2000 = 
(54.625042°, -35.451722°). We adopted the NGC 1399 
centre coordinates from HyperLed43: (a, <5)j20oo = 
(54.621208°, -35.450667°). 

Then we used i ts parametrized i?-band light profile from 
iDirsch et al.l |2003t ) and projected distances of all our UCDs 
to scale the NGC 1399 spectrum, which was then subtracted 
from UCD spectra. The NGC 1399 spectrum was smoothed 
using a fa-spline with equidistant nodes every 4 A. This sup- 
presses additional noise while subtracting it from UCD spec- 
tra and, on the other hand, keeps enough spectral resolution 
not to smear absorption lines in the galaxy spectrum. We es- 
timate the contamination of UCD spectra by comparing the 
median fluxes in the scaled NGC 1399 spectrum and origi- 
nal UCD data. The results (estimated fractions in per cent 
from the total flux) are provided in Table [2] The strongest 
contaminated objects are F-13, F-17, and F-ll. 

We deduced kinematical and stellar population parame- 
ters by fitting the combined FLAMES-HR9 spectra of UCDs 
corrected by the contamination of the NGC 1399 halo us- 
ing the NBursts technique with a grid of simple stellar 
popul ation (SSP) models computed using the lKroupa et all 
(| 19931 ) stellar initial mass function (IMF). In addition to 
the best-fitting values of radial velocity, velocity dispersion, 
SSP-equivalent age and metallicity which are obtained in a 
single minimization loop, for every UCD we computed a con- 
fidence m ap in the age-metallicity sp ace similar to those pre- 
sented in (|Chilingarian et al.ll2008bl ). We notice, that when 
the NBursts technique is used in the x 2 - m appmg mode, i.e. 
fitting a single SSP with fixed age and metallicity at every 
grid node, the algorithm "degenerat es" into the form totally 
equiva lent to the ppxf procedure bv lCappellari fc Emselleml 
l|2004f ) with a single template spectrum. Then we interpret 
the obtained x 2 map in terms of stellar populations. In Ap- 
pendix A we present the numerical experiment aimed at 
quantifying the spectral information contained in different 
absorption features present in the spectral range of our ob- 
servations. We show that the full spectral fitting technique 
in our spectral range using high-resolution pegase.hr mod- 
els successfully breaks the age-metallicity degeneracy and 
can be used to estimate stellar population parameters of old 
stellar populations from FLAMES /Giraffe HR09 data. 
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T able 1. Coordinates of I 
of lMieske etaD l|2008l l. 

id id* 



UCDs excluded from the final sample 



R.A.(J2000) Dec.(J2000) 



F-13 
F-20 
F-28 
F-31 
F-46 
F-62 



FCOS 2-2106 

FCOS 0-2032 
gc319.1 



03:38:29.16 
03:38:56.21 
03:38:25.05 
03:38:19.78 
03:38:30.22 
03:38:49.85 



-35:27:19.9 
-35:24:48.9 
-35:29:25.2 
-35:23:39.5 
-35:21:31.3 
-35:23:36.0 



* ide ntification from 
from iBergond et al.l 



Mieske et alJ l|2004h for F-28 and F-46 and 
20071 ) for F 62. 
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Figure 2. Two-component stellar population fitting of the F- 
13 spectrum. The observed spectrum is shown together with its 
best-fitting template and its decomposition into two components, 
NGC 1399 (green) and a UCD (red). 



In principle, UCDs may contain a mix of different stellar 
populations so that their spectra will not be well represented 
by SSP models. The NBuRSTS technique can be used to fit 
multi-component stellar population models, however, quite 
low signal-to-noise ratios of our spectra did not allow us to 
do it for every object. It was indeed possible for a few bright- 
est ones, such as F-19, were we attempted to fit two SSPs. In 
all cases, the light fractions of young (or intermediate) pop- 
ulation were found to be zero, while old components were 
identical to the single SSP cases. 

Our final sample i ncludes 19 out of 23 UCDs presented 
in iMieske et al.l |2008l ). Two of the four excluded objects, 
F-12 and F-59 had obvious artefacts in the sky subtraction 
due to their very close position to SimCal fibers, while F-60 
and F-64 had too low signal-to-noise ratios in addition to 
difficulties with the sky subtraction. We are adding 6 objects 
to the sample of lMieske et al.l (|200&t ) for which the data were 
obtained, but they were dropped from the final selection, F- 
13, F-20, F-28, F-31, F-46, and F-62. Their coordinates are 
given in Table Q] Four of them, F-20, F-28, F-31, and F- 
46, have low metallicities (see next section). None of them 
but F-13 has archived HST data available. For F-13 it is 
impossible to perform surface photometry with reasonable 
precision, because it is located 0.3 arcmin in projection from 
the photometric centre of NGC 1399 and therefore it is pro- 
jected on the bright part of its spheroid resulting in heavy 
light contamination. 

The approach we use to correct UCD spectra for the 
NGC 1399 halo contamination may introduce biases in case 
of: (1) highly eccentric or distorted isophotes of NGC 1399; 
(2) strong metallicity and/or veloci ty dispersion gradie nts in 
the NGC 1399 halo. According to iDirsch etafl l|2003T ). the 



structure of the NGC 1399 halo is very regular; the veloc- 
ity dispersion profile of NGC 1399 flattens out at a l evel of 
about 230 km s _1 beyond 15 arcsec llSaglia et al.ll2000l) drop - 
ping down only beyond 500 arcsec 1 Schuberth et al.l 120101 ) 
where its surface brightness is too low to cause any con- 
tamination to the CSS spectra . The optical colour profile is 
quite steep |Dirsch et al.l [20031 ) which may be indicative of 
a metallicity gradient. However, the effects of age and dust 
may also play a role. 

In order to test whether our halo correction results are 
strongly affected by radial gradients in NGC 1399, we fitted 
the uncorrected spectra of the three heavily contaminated 
UCDs mentioned above with the model including two SSPs 
with different line-of-sight velocity distribution correspond- 
ing to a UCD and the NGC 1399 halo. At first, we left free 
all kinematical and stellar population parameters of the two 
components. The fitting results for F-13 are shown in Fig [2] 
Then we repeated the test fixing the velocity dispersion of 
NGC 1399 to 230 km s" 1 . In both experiments, the UCD 
kinematical and stellar population parameters, as well as the 
recovered UCD spectra were in excellent agreement. For F- 
17 and F-ll they also precisely matched the values obtained 
from the fitting of halo corrected spectra. The recovered rela- 
tive mass contributions of two SSP models obtained from the 
two-component fitting converted into light contributions are 
also in a good agreement with the values estimated from the 
light profile of NGC 1399. However, for F-13, the NGC 1399 
light contribution to the total spectrum recovered from the 
two-component fitting (73 per cent) turns to be lower than 
that determined from the brightness profile (83 per cent). 
This results in a 0.3 dex lower metallicity (-0.2±0.09 dex) 
and 6 km s" 1 higher velocity dispersion (38±2 km s~ x ) for 
the two-component fitting. 

From this test we conclude that for UCD spectra weakly 
and moderately (up-to 35 per cent) contaminated by the 
NGC 1399 halo, the correction technique based on re-scaling 
the NGC 1399 spectrum according to its brightness profile, is 
consistent with the direct two-component fitting of spectra. 
However, the latter approach becomes rather unstable for 
weakly contaminated spectra because the low overall signal- 
to-noise ratio of our data can not be compensated by high 
spectral resolution for the NGC 1399 component due to its 
very high intrinsic velocity dispersion. 



2.2 Imaging data 

Using Virtua l Observatory t ools, we identified 8 UCDs listed 
in table 3 of lMieske et alJ (|2008l ) as objects without avail- 
able HST photometry. For two of them, F-2 and F-8, F606W 
band archival HST WFPC2 F606W band images available 
through the Hubble Legacy Archive (HLA jf| were found and 
analysed. The data discovery and access procedure com- 
prised several steps: 



(i) The reduced and combined 



FLAMES 



dataset was 

loaded into the VO-Paris Euro3D client^, a dedicated 
Virtual Observatory tool to work with multi-object and 3D 



3 http://hla.stsci.edu/ 

4 http:/ /voplus. obspm.fr/~chil/Euro3D/ 
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Figure 3. FLAMES/Giraffc HR9 spectra, their best-fitting templates (Kroupa IMF), fitting residuals and confidence levels of the age 
and metallicity determinations (inner panels). Spectra, best-fitting templates, and fitting residuals were smoothed with a box of 7 pix 
for clarity. The flux uncertainties were reduced by a factor of 2.65 correspondingly. 



spectorsopic datasets (|Chilingarian et all l2008af ) . The cur- 
rent version of the VO-Paris Euro3D client includes the 
native support of the FLAMES FITS data format. 

(ii) Then, the fibre coordinates were sent to CDS ALADIN 
()BonnareI et alJlifJOoh from the VO-Paris Euro3D client 
using the PLASTIC0 application messaging protocol and then 
displayed on-top of a DSS2 image of the central part of the 
Fornax clu ster. The 8 objects reported not to have HST 
imaging by iMieske et all (|2008l ) were selected from the list. 

(iii) Then, a footprint, i.e. a union of all fields of views of 
all HST imaging observations in a given area of the sky were 
requested from the HLA and displayed using CDS ALADIN. 
Two out of eight objects were highlighted as those covered 
by the footprint. 

(iv) Selecting those two objects one by one in the dis- 
play of CDS ALADIN caused the entries corresponding to 
the WFPC2 field containing them to highlight in the HLA 
Simple Imaging Access query response metadata tree, i.e. a 
structured list of HLA datasets. 

(v) Then, the corresponding WFPC2 images were re- 



trieve d from the HLA and analysed using the GALFIT soft- 
ware (|Peng et al.ll2002j ). 



We used the TinyTim softwar^f] to generate the 
WFPC2 PSF at the positions of F-2 and F-8 in the cor- 
responding WFPC2 fields, which was then used during the 
light distribution fitting to convolve the input model light 
profiles. 

Bo th F- 2 and F-8 were fitted well by single-component 

ISersid 1 19681 ) profiles. For F-2 we obtained the follow- 
ing parameters in the F606VF photometric band: mtot = 
19.98 mag, R e = 0.15 ± 0.01 arcsec = 14 ± 1 pc, n = 
4.9 ± 0.1. For F-8, the parameters are: mtot = 20.31 mag, 
R e = 0.07 ± 0.01 arcsec = 6 ± 1 pc, n = 1.2 ± 0.1. 



5 http:/ /plastic. sourccforge.net/ 



6 http: / /www. stsci.edu/software/tinytim/ 
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Figure 4. Continuation of Fig [3] 



3 RESULTS AND DISCUSSION 

3.1 Internal dynamics and stellar populations. 

The results of the full spectral fitting are presented in Ta- 
ble [2] We provide radial velocities, internal velocity disper- 
sions, SSP-equivalent ages and metallicities for all 24 UCDs 
in our sample. We obtained the mass-to-light ratios of stel- 
lar populations from PEGASE.2 SSP models corresponding 
to the same ages and metallicities. We limited the grid of 
the stellar population models to the maximal age of 15 Gyr. 
There is a number of UCDs where the best-fitting age values 
corrspond to this limit. Although this value slightly exceeds 
the presently adopted age of the Universe, one has to keep 
in mind that there is a well-known "zero-point" problem 
in the stellar population models related to our insufficient 
knowledge of stellar evolution (see "the white paper" by 
IWorthevll2009l for discussion). Therefore, published ages of 
early-type galaxies formally estimated using stellar models 
often slightly exceed the Hubble time. 

We do not compare our stellar M/L ratios with those 
obtained using different stellar population models available 
in literature because in order to get consistent results, one 
has to use the sets of models based on the same ingredients 
(e.g. stellar evolutionary tracks, IMFs, binary fractions) to 



estimate stellar population parameters from the spectra and 
M/L ratios from these parameters. Most ingredients are the 
same for the pegase.hr and pegase.2 models. Even though 
stellar libraries differ, the M/L ratios of the two sets match 
in the U-band. At the sam e time, using different models (e.g. 
iBruzual fc Cha riot 2003) for the M/L ratio estimates would 
require to use them as well for the determination of stellar 
population parameters with the full spectral fitting. How- 
ever, fitting low resolution SSP models against our data in 
a narrow spectral range of the FLAMES HR09 setup would 
not allow us to make any sensible estimates of age and metal- 
licity. 

The uncertainties of the stellar population parameters 
range from 0.03 dex in [Fe/H] and 15 per cent in age for F-19 
(UCD 3) to >0.3 dex in metallicity and totally uncertain age 
for the faintest representatives of our sample such as F-62. 
Generally, there is a clear correlation between the quality 
of age determination and the overall metallicity of a galaxy, 
which can be easily explained as in more metal-rich objects 
absorption lines are stronger, thus better constraining the 
fitting procedure for a given mean signal-to-noise ratio in 
the continuum. 

Thanks to the high spectral resolution of our data and, 
therefore a well sampled line-of-sight velocity distribution in 
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Figure 5. Continuation of Fig [4] 



most UCDs, the vel ocity dispersion - m etallicity degeneracy 
(see Section 1.3.1 in Chilingarian 2006, for details) intrinsic 
to the full spectral fitting in the pixel space makes very little 
effect on the obtained measurements of velocity dispersion. 
Hence, even in cases where the metallicity measurements 
have very large uncertainties due to low signal-to-noise ra- 
tios, the velocity dispersions still remain well determined. 

Our velocity dispersion valu es sometimes d iffer s ignif- 
icantly from those obtained by iMieske et al.l |2008l ) (see 
Fig. [7]), therefore their results of dynamical modelling, i.e. 
aperture corrections, central and global velocity dispersions, 
and, consequently, dynamical mass-to- light ratios have to be 
corrected. We repe ated the dynamica l modelling in exact! 
the same way as in lHilker et all (^OPTf ); IMieske et all (|20Q; 
using new velocity disper s ion es timates and structural prop 
erties from IMieske et al.l (|2008l ) for 13 UCDs, and the pho 
tometric properties of F-2 and F-8 presented above. 



3.2 Dark matter content and stellar mass 
function. 

Using dynamical and stellar mass-to-light ratios, we esti- 
mate the dark matter content of 14 UCDs in our sample 
as ((Af/L)dyn,v - (M/L)*v)/(M/L)dyn,v- The correspond- 



ing fractions in per cent are presented in the last column 
of Table [5] For 9 our of 14 objects the derived dark matter 
contents are consistent with zero, while for F-5, F-U, _F-34, 
and f-53 the dynamical masses turn to be lower than the 
stellar ones resulting in "negative" formally computed dark 
matter fractions, i.e. —100 per cent corresponding to the 
stellar M/L ratio being twice as high as the dynamical one. 
What is the reason for derived negative dark matter frac- 
tions? We either underestimate the dynamical mass and/or 
overestimate the stellar one. 

Underestimatation of dynamical masses may originate 
from imperfect surface photometry (too small effective ra- 
dius) or velocity dispersion measurements biased towards 
lower values. Although our UCDs have quite small spatial 
sizes, most of them are sufficiently well resolved to rule out 
the former possibility. The latter option is also hardly pos- 
sible, because our data have sufficiently high spectral res- 
olution to reliably measure velocity dispersions down to 
5 km s _1 . 

Stellar masses may become overestimated if the adopted 
stellar mass-to-light ratio is too high. The only parameter 
which may strongly affect them at our regime (i.e. old ages, 
intermediate and low metallicities) is the shape of the stel- 
lar mass function at low masses, where stars do not con- 
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Table 2. Projected distances, estimated contamination of UCD spectra by the NGC 1399 light (Ch), Internal kinematics, stellar 
populations, corrected dynamical and stellar mass-to-light ratios and estimated dark matter content of 24 UCDs. 
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Figure 6. Ratio of dynamical to stellar masses vs dynamical 
relaxation time at the half-light radius for 14 UCDs with deter- 
mined dynamical masses. The dotted line denotes the equality of 
dynamical and stellar masses for the assumed stellar mass func- 
tion. Objects above the line either contain dark matter, or have 
either a bottom- or top-heavy IMF; objects below the line must 
have their stellar mass functions affected by dynamical evolution. 



tribute much to the total light, but do change significantly 
the mass because of strongly non-linear behaviour of the 
mass-luminosity relation for main sequence stars. Chang- 
ing the low mass slope power law index from 2.3 to 1.3 (i.e. 
Salpeter to Kroupa) decreases the mass-to- light ratios of old 



stellar populations by ~50 per cent. Does this suggest IMF 
variations in the UCDs of our sample? 

We have to keep in mind that CSSs observed to- 
day might have experienced dynamical evolution ef- 
fects on their stellar mass functions, i.e. the observed 
mass function may differ from the IMF. It is known 
dSpitzerl 1 19871; iBaumgardt fc Makinol 120031 ; iKhalisi et all 
120071 ; iKruiissen fc Mieskd I2009T ) that in globular clusters 
the dynamical evolution causes mass segregation, i.e. mas- 
sive stars are moving toward the centre while low-mass 
stars migrate to the cluster outskirts, where they are tidally 
stripped during the passages close to the centre of a host 
galaxy or through its disc. This creates a deficit of low 
mass stars in a cluster, changing the shape of its inte- 
grated stellar mass function. The characteristic timescale 
of this process is related to the dyn amical relaxation ti me, 
which can be estimated for a CSS (|Mieske et al. I [20081 ') as 



d yn ri/0.0045 Myr, where M dyn is in 



crolax - l og M dy 

Solar masses and r e is in pc. 

In Fig|5]we present the ratio between dynamical (Md y n) 
and stellar (M*) masses versus the relaxation time. We see 
the trend that objects where M* estimates exceed Md yn have 
relaxation times shorter than the Hu bble time. This tren d 
is similar to that presented in Fig 12 of lMieske et al. I (|2008h . 
but in our case the dynamical mass-to-light ratios are nor- 
malised by the stellar ones. It suggests that for some of the 
least massive UCDs dynamical evolution is sufficiently ad- 
vanced to have experienced preferential loss of low mass stars 
and a corresponding flattening of the low-mass stellar mass 
function, such that we overestimate stellar masses assuming 
it to have the Kroupa IMF shape. 

We notice that at the same time very little, if any, dark 
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mass beyond a canonical Kroupa IMF is required to explain 
the dynamical M/L ratios of the more massive UCDs inves- 
tigated here. Since the massive UCDs described in this study 
are dynamically un-evolved, the mass segregation and possi- 
ble tidal stripping of low-mass stars should not affect them. 
T herefore, we can rule out very bottom heavy IMFs like 
in ISalpete 3 ljl955h . which would correspond to 50 per cent 
larger stellar masses and hence imply negative dark matter 
fraction for all investigated sources. 

There have been various claims in the recent litera- 
ture regarding a possible variation or invariance of the IMF 
at cosmological dist ances. While ICappellari et aL I ((Hi); 
iFerreras et al.l (|200ct ) argue against a bottom heavy Salpeter 
type IMF at high redshift, other gravitat ional lens results 
do favour a bottom-heavy Salpeter IM F |Treu et al.l l20ld ; 
lAuger et al.ll2010l ; iGrillo fc Gobatll2010h . Our result for the 
investigated Fornax UCDs is hence in line with the for- 
mer studies. A caveat is that the M/L ratios of UCDs ap- 
pear to vary system atically amongst different environments 
jMieske et allhoOSf l. The case may indeed be different for 
Virgo UCDs which show on average some 40 per cent higher 
M/L at comparable metallicity and age, and which hence 
appear good c andidates for an IMF that deviates from the 
Kroupa form (|Dabringhausen et al.1120081 ). 
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Figure 7. Compa rison of velocity dispersion estimates from 
iMieske et all (I2008T) to our present work. 



3.3 Comparison with literature 

Three objects in our sample, F-l (UCD 2), F-19 (UCD 3), 
and F-24 (UCD 4), have published mea surements of inter- 
nal ki nematics and stellar populations (|Chilingarian et all 
20081)) obtained using exactly the same data analysis tech- 
nique, i.e. the NBuRSTS full spectral fitting FLAMES-LR04 
spectra having similar although slightly wider wavelength 
range but twice lower spectral resolution. The velocity dis- 
persion measurements agree remarkably well for all three ob- 
jects. The metallicity measurements agree within a few hun- 
dredths dex for UCD 3 and UCD 4, however being discrepant 
by ~ 0.4 dex for UCD 2. For t his object, the age estima- 
tions also differ: intermediate in lChilingarian et al.l (|2008bT ) 
and old in our present study. We notice, however, that the 
di screpancy of the UCD 2 age and metallicity measurements 
in lChilingarian et al.l |2008bT ) and our present study follow 
the age-met allicity degeneracy. Given quite poor data qual- 
ity and large size of stella r population confidence levels in 
IChilingarian et all (|2008bl ) we conclude that the discrepancy 
between the measurements can be explained by statistical 
effects. Age determinations for UCD 3 and UCD 4 agree 
between the two studies within 2a of their statistical uncer- 
tainties. 

It is worth mentioning, that the meta l licity estimates for 
all three UCDs presented in lMieske" et all (|2006ft are system- 
atically lower by 0.20. . . 0.25 dex compared to our present 
measurements. 

In Fig [7] (black data points) we present the comparison 
of published velocit y dispersion measurements for 19 UCDs 
(|Mieske et al. I [20081 ) with those obtained in our study with 
the NBuRSTS full-spectral fitting technique. Although the 
general trend agrees, the measurements for individual ob- 
jects are often notably discrepant. 

The reasons for the discrepancy are: (1) template mis- 
match during cross-correlation due to the metallicity differ- 
ence between UCDs and u Cen giant stars served as tem- 



plates; (2) slightly different wavelength ranges used for the 
data analysis (inclusion of the Mg6 triplet region in our 
study); (3) our correction for the co ntamination o f UCD 
spectra by the NGC 1399 halo which IMieske et ail (|2008l ) 
did not apply. 

The metallicity difference between the spectra and tem- 
plates used to analyse them may lead to bi ased velocity 
dispe rsion measurements (see Section 1.3. 1 in IChilingarianl 
l2006f ) at least if the analysis is done in the pixel space, 
thus affecting both NBuRSTS full spectral fitting and FX- 
COR cross-correlation measurements. The low metallicity of 
a template star resulting in shallower absorption lines may 
be compensated by decreasing the velocity dispersion, i.e. 
smearing absorption lines to a lower degree than it should 
be in order to match the line depth in the target spectrum 
being analysed. Then we would expect to see the correlation 
between UCD metallic ities and differences of velocity disper- 
sion measurements in lMieske et alj (|2008l ) and our present 
study, which we do not detect at a statistically significant 
level. This can be explained because the described degener- 
acy between metallicity and velocity dispersion mostly af- 
fects the data for targets with velocity dispersions similar 
to or lower than the instrumental spectral resolution. In our 
case, most low-<r targets have low metallicities well corre- 
sponding to those of the ui Cen template stars. On the other 
hand, massive metal-rich UCDs with relatively high velocity 
dispersions are bright, hence their spectra have good signal- 
to- noise_£atiojM-educing the degeneracy effects. 

IMieske et ail (|2008T ) had to exclude the Mgfe triplet re- 
gion from their analysis as it seemed to bias the velocity dis- 
persion measurements obtained using the cross-correlation 
technique. However, Mgb is the most prominent spectral 
feature in the wavelength range of our spectra, thus con- 
taining a large fraction of spectral information. We per- 
formed the NBuRSTS spectral fitting in the wavelength range 
Arostframc > 5200 A and compared the measurements of ve- 
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locity dispersion with those obtained from the fitting in the 
entire available wavelength range aimed at checking whether 
the full spectral fitting technique also suffers from similar 
biases. The values turned to be consistent within their un- 
certainties, however, being almost half as precise in case of 
the truncated wavelength range. Therefore, we conclude that 
the biases of velocity dispersion measurements obtained by 
cross-correlation including the Mgfc triplet in the wavelength 
range probably originate from the template mismatch when 
using stellar spectra as references, which is minimized in 
our case by selecting the best-matching SSP from the grid 
of stellar population models. As far as these biases affect 
only measurements of low velocity dispersions in objects like 
globular clusters made on relatively high resolution spectra 
and do not seem to show up in studies of relatively massive 
galaxies, we suppose that this mismatch between the spec- 
tra of individual stars and unresolved stellar populations 
originates from subtle absorption-line features and becomes 
important only at high spectral resolution. 

Fin ally, we repea t ed th e velocity dispersion measure- 
ments of lMieske et all (|200ct ) using the FXCOR task in NOAO 
iraf, but now applied to the UCD spectra corrected for 
the contamination of N GC 1399. Using the same stellar 
templates of ui Cen as in lMieske et al.l |2008i ). the discrep- 
ancy to the results obtained with the NBuRSTS technique 
has dramatically decreased which is clearly seen in Fig [7] 
(red data points). We obtain a good agreement between 
the two datasets, with only a small residual systematic off- 
set in the sense that FXCOR dispersions are slightly larger 
than Nburst dispersions for low dispersion values. Thus, 
we conclude that the contamination of the spectra by the 
host galaxy halo is the main reason for biases of the esti- 
mated kinematical parameters, and it has to be taken into 
account in all studies addressing the internal dynamics of 
compact stellar systems. 



3.4 Faber— Jackson and metallicity— luminosity 
relations. 

In Fig. [8] we present the relation between abso- 
lute magnitudes and measured velocity dispersions (top, 
iFaber fc Jacksonlll976l ) and metallicities (bottom) for early- 
type galaxies and CSSs. The giant galaxies are repre- 
sented by 140 ea rly-type galaxies in th e Coma cluster 
from SDSS DR7 (jAbazaiian et al.l 120091 ) obtained using 
the NBurst s technique and vel ocity dispersion measure- 
ments from iBender et all (1 19921). The data for dE/dSOs 



are from I Chilingarian et all (j2008d ) and lChilingarianl (|2009l ) 
(global "main gal actic body" va lues). The local group dwarf 
spheroidals from iMated (| 19981 ) are at the low-luminosity 
end of the diagrams. The data for cE galaxies in cluding 
"legacy" objects are from IChilingarian et all J2009|). tran- 



sition al cE /UCD galax i es are from IChilingarian fc Mamonl 
(|200gl ) an dlPrice etail d2009l ) and t hree additional UCDs 
are from IChilingarian et all J2008bl). The globul a r clus - 
ter data come f r om [McLaughlin fc van der Marell (|2005l ); 
iDiorgovski et all (jl997T ): iMartini fc Hoi (|2004l ) for Milky 
Way, M 31, and NGC 5128 clusters respectively. 

In the upper panel of Fig. [8j we see that all new kine- 
matical measurements of UCDs follow the faint magnitude 
extrapolation of the Faber-Jackson relation (L oc a 4 ) de- 
fined initially for giant early-type galaxies and shown as a 
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Figure 8. Fabcr-Jackson (top) and metallicity-luminosity (bot- 
tom) relations for early type galaxies and compact stellar systems. 
The data sources are given in the text. The dashed line in the up- 
per panel denotes the Faber-Jackson relation (L oc <r 4 ) defined 
for massive early-type galaxies. 



dashed line in the upper panel of Fig.[8j Most cEs and transi- 
tional cE/UCD galaxies reside on it. At the same time, dEs, 
dSOs and the dark matter dominated dSph galaxies deviate 
from the Faber-Jackson relation exhibiting lower velocity 
dispersions for the same luminosity. 

Globular clusters form a different sequence on this di- 
agram having a steeper slope (L oc a 2 ) than the Faber- 
Jackson relation. UCDs of our sample populate the low- 
luminosity extension of the Faber-Jackson sequence down 
to luminosities where it joins with that of globular clusters 
creating a smooth transition. 

On the metallicity-luminosity diagram the situation is 
notably different: the Z — L relation denoted by a dashed 
line in the lower panel of Fig.[HJ is formed by dSph, dE/dSO, 
and giant early-type galaxies, whereas the compact galaxies 
are generally more metal-rich and do not show such a corre- 
lation between their metallicities and luminosities filling the 
entire region above the Z — L relation up-to the solar metal- 
licity. The offset metallicity with respect to the extended 
dwarf spheroidals is consistent with a more efficient self- 
enrichment, and/or them being formed from already pre- 
enriched material in the course of wet galaxy mergers. 
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3.5 Origin of UCD galaxies 

We conclude this paper with a brief discussion on the ori- 
gin of UCDs. One can break down the currently discussed 
formation channels (see Introduction) into two concepts: 

(i) UCDs are tidally stripped remnants of more extended 
galaxies, hence of galaxian origin. As such they would also 
trace the tidal disruption of low-mass dwarf galaxies. 

(ii) UCDs are massive star clusters - or mergers thereof 
- whose formation is closely linked to the formation of the 
bulk of the globular cluster population. 

Regarding the first concept, IChilingarian et al.l (|200g| ) 
demonstrated that a Milky Way sized disc becomes tidally 
stripped on a timescale of 0.5-1 Gyr by an intermediate 
mass cD galaxy similar to the Virgo cluster M 87. Since 
dwarf galaxies usually have lower stellar disc densities, the 
threshing should be even more efficient and quick. A single 
passage near the cD galaxy should be sufficient, suggesting 
a characteristic time of the process to be about 200 Myr. 
This explains why we do not observe many dEN galaxies 
during the process of tidal stripping, this phase is short and 
therefore it is statistically unlikely to observe it. 

The discovery of the first transiti onal cE/UCD galaxy 
M59cO l|Chilingarian fc Mamonll200Sft an d a p opulation of 
galaxies in the Coma cluster ( Price et al.ll2003 ) completely 
filling the mass gap between cEs and bright UCDs, gives 
additional support to the common evolutionary path fol- 
lowed by members of these two classes. However, classical 
cE galaxies can hardly be formed as tidal superclusters in 
galaxy mergers because (1) they have very significant masses 
exceeding 10 9 A/0 which are quite problematic to assemble; 
(2) all nearby cEs rotate and contain either confirm ed cen- 
tral supermassive black holes (|Davidge et al.l [20081 ) or cen- 
tral bumps in th e velocity dispersion distribut ion probably 
caused by them (|Chilingarian fc BergoncfeoiOh . 

Thus, if we accept (massive) UCDs and cEs to share 
the same evolutionary scenario, they should have formed by 
the tidal threshing of more massive progenitors. The clear- 
est candidates for such tidally stripped nuclei are the two 
respecti ve most massive UCDs in the Fornax and Vir go clus- 
ter (e.g. lEvstigneeva et af]|2007l ; Iffilker et all 120071 ). which 
have extended low-surface brightness envelopes. They have 
masses M ~ 1O 8 M0 and clearly stand out among all other 
UCDs in terms of their sizes and luminosities, see also the 
location in the M— a plane in Fig. [8] for the most massive 
Fornax UCD. 

In contrast, the bulk of fainter, lower mass UCDs show a 
smooth luminosity and size transition towards the regime of 
ordinary globular clusters, and their number counts are typ- 
ically well accounted for by an e xtrapolation of the globular 
cluster luminosity function fe.g. iMieske et al]|2004l ). This is 
consistent with a scenario where the bulk of these UCDs 
(in particular the metal-rich ones) formed when the major- 
ity of the red GCs and the spheroids of their host galaxies 
were created (either via a monolithic collapse or multiple 
wet mergers). 
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APPENDIX A: AGE AND METALLICITY 
SENSITIVE SPECTRAL INFORMATION IN 
THE FLAMES/GIRAFFE HR09 SETUP 

Using the technique described in IChilingarianl (|2009l ). we 
quantified the age and metallicity sensitive information in 
the absorption-line spectra at intermediate spectral resolu- 
tion R = 10000 in the wavelength range corresponding to the 
HR09 setup of FLAMES /Giraffe for a single SSP fitting. We 
explored the relative sensitivity of different spectral features 
to the stellar population parameters de fined in the rotated 
age-met allicity coordinate system r\ — 9 (| Chilingarian et all 
l2008d ; IChilingarianl [20091 ) with the 77 axi s parallel to the 
age-met allicity degeneracy |Worthevlll9"94T ): 

r? = (3Z + 21og 10 t)/Vl3; 

0= (-2Z + 31og 10 t)/\/l3 (Al) 

The contribution of every pixel at every wavelength to the 
total x 2 depends on many parameters. In order to estimate 
the fitting procedure sensitivity S(X,po, ■ ■ ■ ,p n ) to a given 
parameter p; at a given wavelength A, we should compute 
the corresponding partial derivative of the template grid re- 
moving the global continuum shape, i.e. accounting for the 
multiplicative continuum variations. Hence, the response to 
a stellar population parameter r\ (and similarly, 9) at a given 
point (to, Zo,<Jo) of the parameter space is expressed as: 

S(X,rj)\(t ,z ,v ) = 

^- x 2 {P lp (\, v , 9)\ {to , Zo) (T(A, t, Z) ® £(a„))}, 

(A2) 

where £ is the LOSVD, T(A, t, Z) is the flux of a given SSP 
template spectrum characterised by its age and metallicity 
at a given wavelength, Pi p is the p-th order multiplicative 
Legendre polynomial. The meaning of this equation is that 
the squared fitting residuals of a spectrum (770 + Ari,9o) 
against the model (rjo,6o) at every pixel would correspond 
to its contribution to the overall \ 2 when varying r\. 

In practice, we compute these derivatives numerically 
using a "one-side" approach, i.e. computing the function val- 
ues at a point of interest and at a nearby point slightly off- 
set from it on a given coordinate. For our test we chose two 
SSP models with the age t = 12 Gyr and metallicities —1.0 
and —0.4 dex correspondingly. Then we convolved them 
with Gaussian kernels corresponding to the internal veloc- 
ity dispersions of 10 and 20 km s -1 and, thus, obtained four 
model spectra or "reference SSPs". Then for every of them 
we varied the r\ and 9 parameters corresponding to their 
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Figure Al. Spectral distribution of the stellar population sen- 
sitive information for four stellar populations in 5 A-wide bins 
in the FLAM ES/Giraffc HR09 setup. Corresponding age, metal- 
licity, and velocity dispersion are indicated in every panel. The 
studied SSPs arc shown in black; blue and red histograms display 
the relative importance of spectral bins (in per cent) for the de- 
termination of n and 9 stellar population properties respectively. 
Green bars denote the bands defining the following Lick indices: 
Mgb, Fe5270, and Fes335 with yellow side-bars corresponding to 
the corresponding pseudo-continuum definition regions. 



ages and metallicities (see Eq. lAl[l by 0.05 dex, which corre- 
sponds to (12.8 Gyr, -0.96 dex) and (13.2 Gyr, -1.03 dex) 
for the metal-poor models and to (12.8 Gyr, —0.36 dex) 
and (13.2 Gyr, —0.43 dex) for the metal-rich ones. Later, 
we fitted these models against their "reference SSPs " us- 
ing the ppxf procedure (|Cappellari &i Emsel lcm 2004) with 
the 10th order multiplicative polynomial continuum. The 
fitting residuals obtained by this procedure correspond to 
the quantity defined in Eq. IA2I 

We have co-added the information on r/-, and 6- 
sensitivity in 5 A bins in the wavelength range between 5100 
and 5380 A and normalised it by the total value of non- 
reduced x 2 thus obtaining the relative importance of every 
5 A-wide bin to the determination of the stellar population 
parameters. The results for our four "reference SSPs" are 



shown in Fig I All One should not directly compare blue and 
red curves, because they represent the normalised quanti- 
ties. The absolute values of S(X, rf) as defined by Eq. IA2I 
are 2-2.6 times higher than S(X, 9) explaining the elongated 
shapes of the 1-rj uncertainty ellipses in the age-metallicity 
space. 

The age and metallicity response cannot be determined 
directly using this approach, because these two parameters 
are degenerated. Then, the formally computed sensitivity 
will become overestimated. 

We see that ^-sensitive information (the one allowing us 
to break the age-metallicity degeneracy) and f?-sensitive in- 
formation represented by red and blue histograms in Fig. lAll 
respectively are distributed slightly differently along the 
wavelength range. In general, the blue and red histograms 
are much more similar to each other than i n the case of 
lower spectral resolution as presented in fig. 2 of lChilingarianl 
(|2009l ). However, we see in Fig. lAll that at high spectral res- 
olution the sensitivity across the age-metallicity degener- 
acy (blue histogram) is stronger correlated with equivalent 
widths of metal absorption lines than that along the de- 
generacy (red histogram) which is more bound to the Mgfe 
triplet and fainter absorptions around 5140 A and 5230 A 
but not to the strong iron multiplet defining the Fe527o Lick 
index. 

The main conclusion we draw from this test is that when 
we use the full spectral fitting at a sufficiently high spec- 
tral resolution, then even in the narrow wavelength range of 
the FLAMES '/Giraffe HR09 setup including no "traditional" 
age indicators such as Hf3, we are able to constrain stellar 
population age with a high level of accuracy. 

The whole idea of the stellar population determination 
using the full spectral fitting is related to the difference of 
age and metallicity response functions of different pixels in 
a spectrum. These are related to the flux response of a stel- 
lar spectrum at a given wavelength to effective temperature, 
surface gravity, and metallicity, because we use the stellar 
spectra as principal ingredients to build stellar population 
models. Even if we deal with the spectral region containing 
mostly lines of metals, which become deeper when the stel- 
lar population metallicity is increasing or the population is 
becoming older, the exact behaviour of individual lines (or 
groups of lines) will be slightly different which is again con- 
nected to the absorption line properties in stellar spectra. 
For example, the behaviour of the Mg6 triplet equivalent 
width as a function of T e g- is not the same as that of the 
Fe527o multiplet, while the metallicity dependence is rather 
similar. These differences for different sometimes very faint 
spectral lines and molecular bands are co-added at every 
pixel included in the fitting and at the end bring enough 
contribution to the \ 2 m order to determine the age of the 
stellar population from the fitting procedure. However, if the 
intrinsic velocity dispersion of the stellar system being stud- 
ied is high (or, alternatively, the spectral resolution is low), 
then the age and metallicity sensitive information in faint 
absorption line features is rapidly decreasing, hence increas- 
ing the relative importance of the most prominent spectral 
features as those used to define the Lick system. Therefore, 
while we can successfully determine stellar population prop- 
erties of UCDs from our FLAMES/Giraffe HR09 data, it 
would be virtually impossible to do so for giant early-type 
galaxies, because of their high intrinsic velocity dispersions. 



